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Electronic structure calculations were performed on the compound C0S2 , an itinerant ferromagnet 
whose magnetic properties can be understood in terms of spin fluctuation theory. We have identified 
nesting features in the Fermi surface of the compound, active for long wavelength spin fluctuations. 
The electronic structure of the material is close to a half-metal. We show the importance of in- 
troducing spin-orbit coupling (SOC) in the calculations, that partially destroys the half-metallicity 
of the material. By means of transport properties calculations, we have quantified the influence of 
SOC in the conductivity at room temperature, with an important decrease comparing to the GGA 
alone conductivity. SOC also helps to understand the negative of the material, whose conductivity 
varies by a few percent with the introduction of small perturbations in the states around the Fermi 
level. 

PACS numbers: 



I. INTRODUCTION 

Half metals^ are materials where metallic con- 
duction occurs only through one of the spin chan- 
nels, leading to a perfectly polarized electron cur- 
rent. One of the spin channels presents a gap at 
the Fermi Energy (Ep) but the other one has bands 
crossing it. This makes them an interesting family of 
materials, in particular for applications in the field 
of spintronics^ where spin injection is required for 
controlling the charge and spin currents separately. 
Moving a strongly spin-polarized current is highly 
sought for, and half-metals are a good candidate for 
this. 

Several electronic structure calculations based on 
density-functional theory (DFT) have been done 
in the last few years, describing C0S2 as a half- 
metallic ferromagnet^— However, the electronic 
and magnetic properties of this material are not 
fully understood. Point-contact Andreev reflection 
measurements^ showed a relatively low-spin polar- 
ization of 56% at 4.2K This picture was confirmed 
in the studies by Wang et alH in the Co\- x Fe x Si 
series, where they show that the spin polariza- 
tions can be tuned in the range —56% < P < 
+85%. Brown et alJ£ found from polarized neutron- 
diffraction measurements that half metallicity does 
not occur in the ferromagnetic phase. A sharp pho- 
toemission peak at the Fermi energy in the ferromag- 
netic phase originates in the bottom of e g i subband 
due to the exchange 3, indicating that the e g i band 
is partially filled in the ferromagnetic phased de- 
stroying the half metallicity. 

Recently it was found that the resistivity, spe- 



cific heat and magnetic susceptibility of C0S2 
is dominated by exchange-enhanced spin-density 
fluctuations.— This previous work also show that the 
GGA method^ is enough to describe the electronic 
and magnetic of C0S2 due to the non localized 
nature of the compound. There is no need to intro- 
duce strong correlation effects to describe accurately 
its electronic structure properties but not enough to 
describe the enhancement of, e.g. the value 7 ob- 
tained from specific heat , which is enhanced due to 
spin fluctuations unaccounted for in DFT calcula- 
tions. 

The important role of spin-orbit coupling in a re- 
alistic description of half-metallic materials has been 
discussed in literature about half-metalsiil Usually 
the coupling between spin-up and spin-down states 
caused by a sizable spin-orbit coupling introduces 
non-negligible minority states at the Fermi level, de- 
stroying the half-metallicity. In general, this will 
depend on the strength of the interaction (that in- 
creases with atomic size) and also on the size of the 
gap and the position of the band edges in the insu- 
lating spin channel. The importance of minimizing 
the effects of spin-orbit coupling on destroying half- 
metallicity when designing new half-metals has been 
discussed in the past,— where utilizing cations with 
completely filled shells was suggested in order to re- 
duce the spin-orbit effects by minimizing the orbital 
angular momenta and, in the case of studying pos- 
sible compensated half-metals, increasing the likeli- 
hood of the materials having no net magnetization. 

However, most works based their studies solely on 
the changes in the density of states at the Fermi 
level, but the degree of half-metallicity will be infm- 



enced by the changes in the conduction properties 
that derive from the reordering of states around the 
Fermi level introduced by spin-orbit coupling. More- 
over, the conductivity of the minority spin chan- 
nel will have an activated component that needs 
to be quantified beyond the value of the density of 
states at . In this paper, we will make an analy- 
sis of the electronic and magnetic properties of the 
ground state of C0S2 and introduce spin orbit cou- 
pling (SOC) to determine its importance in the half- 
metallicity of the material. We will also focus on ex- 
plaining, from a band structure point of view, what 
features help us understand its behavior as an itiner- 
ant ferromagnet, governed by spin fluctuations. For 
analyzing the conduction properties of the material, 
its spin-dependence and the importance of spin-orbit 
coupling, we have made several transport properties 
calculations and have compared them with experi- 
mental measurements. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations were performed 
within density functional theory 19 using wien2k 
software ) 20 i 21 which utilizes an augmented plane 
wave plus local orbitals (APW+lo) 22 method to 
solve the Kohn-Sham equations. This method uses 
an all-electron, full-potential scheme that makes no 
shape approximation to the potential or the electron 
density. The exchange-correlation potential utilized 
was the Perdew, Burke and Ernzerhof (PBE) version 
of the general gradient (GGA). 16 The geometry op- 
timization was carried out minimizing the forces in 
the atoms and the total energy of the system. The 
parameters of our calculations depend on the type 
of calculation but for any of them we converged with 
respect to the A;-mesh and to RmtKmax- Values used 
of the fc-mesh are 7x7x7 sampling of the full Bril- 
louin zone for electronic structure calculations and 
geometry optimization, 40 x 40 x 40 for the band 
structure calculations and the very fine details of the 
influence of spin orbit coupling in the half metallic- 
ity of the compound. R m tK max = 6.0 is chosen for 
all the calculations. Local orbitals were added for a 
bigger flexibility in dealing with the semicore states. 
Muffin tin radii chosen were the following: 2.31 a.u. 
for Co and 2.04 a.u. for S. All the calculations were 
converged with respect to all the parameters used, to 
the precision necessary to support our calculations 
(converged forces, total energy differences, etc.). For 
the calculations of transport properties we utilize 
the BoltzTraP code^ that uses the energy bands 
obtained from the wiEN2k software. We used a k- 
mesh in a 43 x 43 x 43 sampling of the full Brillouin 
zone, where convergence was achieved. The calcula- 



tion of transport properties requires a fine mesh to 
carry out the integrations in the Brillouin zone. Spin 
orbit coupling (SOC) was introduced in a second- 
variational procedure* 2 ^ 

III. STRUCTURE 

C0S2 crystallizes in a cubic pyrite structure^ as 
shown schematically in Fig. [TJ Co and S atoms 
are located respectively at the Wyckoff positions 4a 
(0,0,0) and 8c (u,u,u) of the space group Pa3, no. 
205. Co atoms are octahedrally coordinated by 6 S 
neighbors. 

From X-ray diffraction experiments on large single 
crystals (3x3x3 mm 3 )^ we obtained the lattice 
parameter of the material to be a= 5.518 A. Also, 
electronic structure calculations were performed to 
optimize ab initio the structural free parameter 
the S position. This was calculated in order to min- 
imize interatomic forces and the total energy of the 
system. The value obtained was u= 0.387 (see Table 
HI, in good agreement with previous works.™ Calcu- 
lations to obtain the lattice parameter ab initio were 
performed yielding a lattice parameter a = 5.55 A, 
which is within the typical accuracy of density func- 
tional theory calculations (2% in volume). 

TABLE I: Crystallographic positions at T = that result 
from our structure optimization. 





Crystallographic 




Atom 


position 


Coordinates 


Co 


4a 


(0,0,0) 


S 


8c 


(0.387,0.387,0.387) 



IV. RESULTS 

A. Electronic structure and magnetism 

The electronic structure of C0S2 is characterized 
by a large covalency due to the extended S p bands 
overlapping largely the Co d bands close to the Fermi 
level. Also, the itineracy of the compound helps in 
making difficult to establish an ionic picture for the 
material. Some earlier reports suggest an effective 
Co 2+ :d 7 stated distinct from the purely ionic pic- 
ture. 

Our results show that the orbital occupations are 
close to the Co 2+ :d 7 state. As expected, a large co- 
valency can be seen in the density of states (DOS) 
plots in Fig. [2] Big spectral weight from S p states 
appears near Ep and also significant weight can be 
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FIG. 1: (Color online) Cubic pyrite structure of the C0S2 
compound. Large blue spheres represent Co atoms and 
small yellow spheres represent S atoms. 
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FIG. 2: (Color online) Spin-up (positive values) and 
spin-down (negative values) total density of states (DOS) 
plots for the C0S2 in the ground state ferromagnetic con- 
figuration. Partial density of states have been plotted 
also (dash blue lines for Co total DOS and dot green 
lines for S total DOS). Fermi energy is at zero. The 
system is half-metallic. 



seen at the top of the conduction bands, showing 
the S p bands are not fully occupied. In this DOS 
plot, calculated within the GGA approach, we see 
C0S2 presents a half-metallic character, as noted in 
several previous worksiii^ It shows the fully occu- 
pied tig bands of about 0.7 eV bandwidth at about 
-1 to -2 eV, that are spin-split by about 0.6 eV, and 
the broad ej bands crossing the Fermi level in the 
majority channel, being 2 eV wide. The e^ bands 
are unoccupied, with its edge lying very close to the 
Fermi level. This proximity will make their occu- 
pancy non-zero under perturbations such as spin- 
orbit coupling, temperature or the application of an 




FIG. 3: (Color online) Charge difference of the major- 
ity and minority states of C0S2 compound. Large blue 
spheres represent Co atoms and orange small spheres 
represent S atoms. The charge difference plotted cor- 
respond to an isosurface at 0.05 e/A 3 produced using 
XCrySDen.— The shape of the spin density isosurface 
represents the orbital that carries the spin moment, be- 
ing a mixture of e 3 states. 



external magnetic held. A small amount of electron 
doping would also be able to populate them. 

The value of the total moment obtained by ab ini- 
tio calculations is 1.00 ^_b/Co in the ground state 
ferromagnetic solution, in agreement with the value 
of the saturation magnetization found experimen- 
tally, of about 0.9 /ib/Co.— The moment comes 
from the partially filled majority e g band. To il- 
lustrate this description of the electronic structure, 
we have produced a differential electron spin density 
plot (Fig. [3]). This shows the difference in electronic 
density between the majority and minority states of 
the compound (yellow), in order to see the spatial 
distribution of the only unpaired electron (charge 
difference plot will show basically the density that 
corresponds to the electron, because the t 2g bands 
are completely full). Charge accumulation can be 
seen along the Co-S direction (like d3 Z 2_ r 2 or d x 2_ y 2 
in an octahedral environment) . The large bandwidth 
of the e g bands (about 2 eV) , much higher than the 
crystal field splitting between the e g sub-levels, leads 
to a similar occupation of both e g orbitals, as we see 
in Fig. El 

Our results show that a calculation at the GGA 
level is enough to capture the features of the elec- 
tronic structure of this itinerant ferromagnet, as 
shown in a previous work.— The introduction of 
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FIG. 4: (Color online) Band structure plot of the: a) majority spin states, b) minority spin states and c) calculation 
including spin orbit coupling with the magnetization along 001 direction. Observe that the curvature of the bands 
remains unaltered when spin-orbit coupling is turned on. Fermi level is at zero 



strong correlation effects by the LSDA+U method 28 
does not provide more information about the elec- 
tronic structure of this itinerant system, in agree- 
ment with the work of Kwon et alJ± The correct 
ferromagnetic ground state is obtained at the GGA 
level (a FM solution is more stable by 75 meV/Co 
than the paramagnetic phase). A peak of the DOS 
at the Fermi energy (not shown) is found in a non- 
magnetic calculation, suggestive of the system being 
close to a ferromagnetic instability according to the 
Stoner criterion 2 ^ (N(Ep)I > 1 to favor ferromag- 
netism). 

Band structure plots (Fig. and b for the major- 
ity and minority spin channel, respectively) can be 
used to understand the essential features of the elec- 
tronic structure of C0S2 . Being a half- metal, all the 
spectral weight around the Fermi level comes from 
the majority spin channel (Fig. H£i). The minority 
spin channel (Fig. Hp) shows a gap of about 1 eV 
at the Fermi level. A big electron pocket, centered 
at the r point, and a large hole pocket appear in 
the T-R direction. Additional electron pockets ap- 
pear near the M point. These results are similar to 
other band structure plots published using the GGA 
scheme^ except for the electron pocket centered at 
the r point, which is absent in their calculations. 

We have also performed calculations including 
SOC, with the magnetization along the main di- 
rections of the cubic structure: (001), (110) and 
(111). Calculations yield an orbital magnetic mo- 
ment 0.037 fiB (quenched by the crystal field^S), sim- 
ilar to that reported by Antonov et al£^ Together 
with the local spin magnetic moment of 0.92/1 ^ (in- 
side the Co muffin tin sphere) we obtain Mi/M s = 
0.04, i.e., < l z > I < s z >= 0.08, in good agreement 



with experimental measurements (O.l)^ 

A very small energy difference is obtained be- 
tween the calculations with the magnetization along 
the three main directions of the crystal, suggesting 
a very small magnetocrystalline anisotropy in this 
largely isotropic material. The energy difference be- 
tween the easy and hard axis is smaller than 50 
/xeV/Co atom. 

Analyzing the band structure (Fig. 0J:) obtained 
from a GGA+SOC (001) calculation, we observe 
that the introduction of SOC leads to small shifts 
in the bands, but the positions of the pockets in 
the Brillouin zone and the curvature of the bands 
remain roughly unchanged (we see all the electron 
and hole pockets discussed above). A splitting in 
the bands corresponding to the electron pocket in 
the T-R direction, and a variation in the size of the 
other pockets can be seen. Below we discuss further 
the quantitative effects in the electronic structure of 
introducing SOC. 

For a better description of the states near the 
Fermi level, Fig. [5] shows the Fermi surface of the 
compound, within a GGA+SOC-based calculation 
with the magnetization set along the (001) direc- 
tion. In the plot we observe a large isotropic elec- 
tron pocket around the T point. Experimentally^ 
the spin susceptibility of this compound can be un- 
derstood in terms of spin fluctuation theory^ For 
this type of systems, peaks in the spin susceptibility 
can be associated to nesting features in the Fermi 
surface. In the case of C0S2, some nesting features 
can be observed in the T — M — R plane (Fig. [5] 
b)), with a nesting vector that we have denoted as 
Q n . The magnitude of Q n would correspond to long- 
wavelength spin fluctuations (A = 6a, with a being 
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FIG. 5: (Color online) Fermi surface of the C0S2 ob- 
tained for the majority spin channel in a FM solution, 
within a GGA+SOC-based calculation along the (001) 
direction, a) and b) show the 2-dimensional Fermi sur- 
face in the T-X-M and F-M-R planes, respectively, c) 
shows the 3-dimensional Fermi surface. In the inset b) 
we have denoted the nesting vector as Q n (which cor- 
responds to a wavelength of approximately 6 times the 
lattice parameter) . This would lead to a peak in the spin 
susceptibility associated with that wave-vector, together 
with the ferromagnetic ones. 



the lattice parameter), that are believed to exist in 
weak ferromagnets such as CoS^ and also MnSi^i 
together with many other systems^ The rest of the 
Fermi surface, i.e. the hole pocket along T-R di- 
rection and the small electron pockets in the X-M 
direction near M, is largely incoherent and will lead 
to negligible nesting. 



B. Half-metallicity 

C0S2 has been predicted to be a good source of 
spin polarized electrons;-^ hence, we will try to ad- 
dress this point with our calculations. GGA alone 
leads to half-metallic behavior (Fig. [2]), as we saw 
above. Due to the proximity of the conduction band 
edge to the Fermi level in the minority spin channel, 
SOC destroys the half-metallicity noticeably, as we 
can see in Fig. |6l No minority DOS is seen in the 
GGA calculation, but magnifying the minority DOS 
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FIG. 6: (Color online) Spin-up (positive values) and 
spin-down (negative values) total DOS plots for the C0S2 
in the ferromagnetic configuration with a GGA (solid 
red line) and a GGA+SOC (SOC along (001) direction) 
(dash blue line). Fermi energy is at zero. The sys- 
tem is half-metallic in the GGA and 1 metallic in the 
GGA+SOC 0. The inset shows a zoom of the minor- 
ity DOS near the level of both solutions (GGA and 
GGA+SOC) 



of the calculation including SOC (in the inset), we 
see the DOS is nonzero at the Fermi level for the 
minority spin channel. This has important conse- 
quences for the conduction properties of the mate- 
rial. We can compare the values of the DOS at the 
Fermi level for the two spin channels, being no longer 
infinite but N(E F ) T /N(E F ) i = 540 when SOC is 
considered, giving a polarization of the DOS at the 
Fermi level P = [N t {E F ) - Ni(Ef)]/[N^(E f ) + 
N±(E F )] of 99.6%, very large still. 

For analyzing further the degree of spin polariza- 
tion of the carriers in this system, we have calculated 
the transport properties from our band structure cal- 
culations using a semiclassical 1 based on Boltzmann 
transport theory^ In order to understand the spin 
polarization of the carriers, it is necessary to go be- 
yond the values of the DOS at the Fermi level. One 
possible way is to compare the conductivities of both 
spin channels. The inset of Fig. [7]shows the temper- 
ature dependence of the ratio between the subtrac- 
tion and the sum of the conductivities coming from 
both spin channels, i.e. P CT = {o uv -a dn )l (a up +a dn ), 
defined analogously to that of the DOS. This quan- 
tity ranges from ±1 for a full-polarization to for 
no polarization and will give us the spin polariza- 
tion of the carriers at a given temperature. When 
no SOC is considered and polarization of the DOS 
is infinite, it is still very large in terms of the con- 
ductivities (spin polarization of the carriers), being 
approximately P CT = 0.998 at room temperature. 
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FIG. 7: (Color online) Percentage (Psoc)of the bet- 
ter conduction of GGA versus GGA+SOC calculations. 
This is calculated by subtraction of the conductivities 
coming from GGA and GGA+SOC calculations, nor- 
malized to the sum of both them (solid orange circles) 
(Psoc= {&nosoc — o- soc )/(<t„ osoc + a soc )). GGA alone 
leads better conduction properties, so the GGA+SOC 
implies larger resistivity. The inset show the 2 of the 
conductivities coming from both spin channels over the 
sum of these conductivities within the GGA (solid red 
diamonds). At room temperature, a up contributes in 
a 99.8 % to the total conductivity and <Jdn 0.2 % con- 
tributes it to the total conductivity in the GGA 0. 



One would naively expect that the introduction 
of SOC leads to a substantially larger conductivity 
coming from the minority spin channel compared to 
the case without SOC, because of the non-zero DOS 
at the Fermi level. To gain some insight on that, we 
present in Fig. [7]the temperature dependence of the 
calculated ratio between the subtraction and the ad- 
dition of the conductivities calculated with and with- 
out SOC: P 

SOC - {p'nosoc ^ soc) i {p^nosoc~^~ ^ soc 

). We 

observe that the solution without SOC has a larger 
conductivity than the solution including SOC up to 
room temperature, both values differing by about a 
few percent. Thus, the bigger change in the con- 
ductivity when introducing SOC in the calculations 
is not the additional contribution from the minor- 
ity spin states, but an overall reduction in the con- 
ductivity. Even though a separation in conductivity 
from both spin channels is not possible when cal- 
culating with SOC, the results suggest a very large 
polarization of the carriers will still be present when 
the calculations include SOC. This effect will ac- 
count for a reduction in the half-metallicity of just a 
few percent. However, to understand the experimen- 
tal observations, one needs to consider other effects 
as possible sources against a full spin polarization, 
such as defects, spin excitations or non-quasiparticle 
state o 35 i 36 to be important. 




50 100 150 200 250 



T(K) 

FIG. 8: (Color online) Evolution of the magnetoresis- 
tance with temperature measured comparing the values 
of resistivity at H = 9T and H — 0T applied magnetic 
field according to the formula MR=100(/59t — Pot)/ Pot- 
Below 25K we observed positive MR, and negative MR 
are obtained in the rest of FM phase. 

In a material with such a large polarization of 
carriers at the Fermi level, a negative magnetore- 
sistance is expected to occur in the FM region of 
the phase diagram. Figure [5] shows the measure- 
ments of the magnetoresistace (MR) in an external 
magnetic field of 9 T, (MR=100(p 9T - Pot)/ Pot)- 
These show a negative MR in a temperature range 
within the FM phase (below the Curie temperature 
T c = 122if 15 ) with a value ranging between -3% and 
-9%. At low temperatures, a large positive value of 
the MR is observed and, above the Curie point, in 
the paramagnetic phase, the MR becomes again pos- 
itive, approaching a constant value of about 3% at 
high temperatures. MR does not reach a very high 
value, consistent with the fact that spin polarization 
of the carriers is not complete, as discussed above. 

V. SUMMARY 

We have analyzed the electronic structure and 
magnetic properties of this itinerant weak ferromag- 
net, dominated by spin-density fluctuations. We 
have related these to the nesting features in the 
Fermi surface associated to long wavelength nest- 
ing vectors. Being an itinerant system, the GGA 
scheme can describe its properties and electronic 
structure accurately, not needing to introduce the 
electron correlations with means of LDA+U. GGA 
calculations alone predict C0S2 to be half metal- 
lic, but when introducing the important effects of 
SOC effects, the results of the calculations show 
SOC destroys the half-metallicity, non-zero occupa- 
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tion (band-crossing) of the minority channel at the 
Fermi level appear. Our transport properties cal- 
culations show that changes in the conductivity as- 
sociated to the introduction of SOC in the calcu- 
lations do not imply a significant break-up of the 
half-metallic behavior, but just a few percent. These 
results are consistent with our experimental observa- 
tion of negative magnetoresistance in the material. 
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